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Abstract 

Description of time evolution of the relativistic unstable electromagnètic system 
consisting of Fermi - Dirac particle interacting with electromagnètic field, in the 
framework of the Liouville space extension of quantum mechanics is done. The 
work was carried out on the basis of Prigogine's unified formulation of quantum and 
kinetic dynamics. The eigenvalues problem for the relativistic Hamiltonian of the 
electromagnètic system was solved. The obtained results can be used as the ground 
for the further studies of the observed physical processes such as bremsstrahlung, 
relaxation of excited states of àtoms and atòmic nuclei, partides decay. 
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1 Introduction 



It is known triat the description of physical world on the basis of fundamental 
classical and quantum theories is defined by the laws of the nature as determin- 
istic, time reversible. The time in the usual formulation of dynamics does not 
have the chosen direction and the future, and the past are not distinguished. 
However, it is also noted that the facts given before are in the contradiction 
to our experience, because the world surrounding us has obvious irreversible 
nature. In this world the symmetry in the time is disrupted and the future and 
the past play different roles. Difference between the classical description of the 
nature and those processes in the nature which we observe creates the conflict 
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situation. I. Prigogine noted that the solution of this problem is impossible 
(at least accurately) on the basis of the conventional formulation of quantum 
mechanics. Therefore one should speak about the alternative formulation of 
dynamics that makes it possible to include the irreversibility in a natural way. 
In this connection the studies of the irreversible processes at the microscopic 
level - the microscopic formulation of the irreversibility represents for me the 
special interest. 

The alternative formulation of dynamics found its embodiment in the works 
of Brussels-Austin group that was headed by I. Prigogine for many years. The 
authors of the approach deny the conventional opinion that the irreversibility 
appears only at the macroscopic level, while the microscopic level must be 
described by the laws, reversed in the time. The mechanism of the asymmetry 
of processes in the time, which made it possible to accomplish a passage from 
the reversible dynamics to the irreversible time evolution one was developed. 
In the approach the law of the increase of entropy is accepted as the funda- 
mental that determines the "arrow of time", the difference between the past 
and the future. Thus, new irreversible dynamics with the disrupted symme- 
try in the time was formulated. In the approach of Brussels-Austin group the 
irreversibility is presented as the property of material itself and is not defined 
by the active role of the observer. From the other side the approach allows to 
solve the problems, which could not be solved within the framework of classi- 
cal and quantum mechanics. For example, now we can realize the program of 
Heisenberg - to solve the eigenvalues problem for the Poincare's non-integrable 
systems, which could not be solved within the framework of traditional meth- 
ods. The approach solves the bàsic problem, designated by Boltzmann and 
Planck - to formulate the second law of thermodynamics at the microscopic 
level. 

The general formalism of Brussels-Austin group approach was developed in 
works [1] - [22]. In these papers the bàsic ideas of the transition from de- 
terministic dynamics to the irreversible description are formulated. In mono- 
graphs [1], [2] the universal survey of ideas and principies of the alternative 
formulation of dynamics is given. A unified formulation of dynamics and ther- 
modynamics is done, for example, in works [4], [5]. It is carried out the study 
of the non-integrable systems, where the alternative formulation of quantum 
mechanics for non-integrable systems is proposed [6], [10]. The " subdynamics" 
approach is developed in works [3], [7], [8], [11]. In paper [15] the method was 
adapted for the explicit computation of the eigenvalues problem for the Renyi 
maps, baker's transformations, Fridrichs model. The eigenvalues problem for 
the Liouville operator L is solved in the framework of a complex, irreducible 
spectral representat ion [17], [18]. The role of the thermodynamic limit for 
Large Poincare systems is investigated in work [22]. In works [23], [24], scat- 
tering theory in superspace and the three-body scattering theory for finite 
times is developed. In [25] - [27] in the framework of Friedrichs model the 
problem of description of quantum unstable states including their dressing is 
investigated. The problem of the complex spectral representation of Liouville 
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operator in the extended Liouville space outside the Hilbert space is solved. 
It is shown that the dressed unstable state described by a density matrix can 
be expressed in the terms of the Gamow vectors. The Gamow vectors are in- 
vestigated also in works [28] - [32]. The formalism determines the operator of 
microscopic entropy M and also the time operator T [33] . The time operator T 
is constructed for a quantum system with unstable particle. In work [34] one- 
dimensional gas with 5-function interaction is examined. Formalism found its 
further development in works [40] - [46] . So in work [40] on basis of Friedrichs 
model a microscopic expression for entropy is obtained. The simple model of 
interaction harmònic oscillator with a field is developed in work [41]. A prob- 
lem of quantum decoherence for a particle coupled with a field was considered 
in work [42]. It was shown that the decoherence in the field is the result of ir- 
reversible process. In work [43] analysis of the Hegerfeldt's theorem is carried 
out. The analysis of the short-time behaviour of the survival probability in 
the framework of the Friedrichs model was done in work [44]. The two models 
of relativistic interaction are examined in the work [45]. The models involve 
two relativistic quantum fields. They are coupled by the simplest cúbic and 
quadratic interaction. A pair of identical two-level àtoms interacting with a 
scalar field are considered in the work [46]. The qüestions of irreversibility are 
developed also in works [47] - [58]. 

At present in the framework of Prigogine's ideas the great number of works 
with the use of different models of interaction was executed. They are the 
Friedrichs model (see [35] - [39]), the models used interaction of simple cúbic, 
or quadratic form as, for example, in works [28], [45] or 5-function interac- 
tion [34]. Therefore, at present moment, it is very interesting and necessary 
to continue further development of the formalism with the use of realistic rel- 
ativistic Hamiltonians. 

In the paper I examine the possibility of application Brussels-Austin group's 
ideas, for the description of time irreversible evolution of the quantum un- 
stable electromagnètic system consisting of Fermi - Dirac particle interacting 
with electromagnètic field. The model of relaxation of the system with photon 
emission is investigated. It is important that the Hamiltonian of interaction 
is determined on the basis of the requirement of the gauge invariance of the 
model. The definition of the interaction model is done in section 2. In sec- 
tion 3 I present the Liouville formalism, "subdynamics" approach. In section 
4 the initial expression for the density matrix is formulated. The task of the 
complex spectral representation of Hamiltonian is solved in section 5. The 
expression for the density matrix describing the evolution of the relativistic, 
unstable electromagnètic system depending on the time is obtained in section 
6. Numerical calculation are given in section 7. 
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2 Definition of the interaction model 



For the operator of fermion field ip(x) we have the following decomposition [59] 



^ <+, W = (^/0' /2 ·' r W^M' i P. (i) 

where c r {p) (cj(p)) is the operator of destruction (creation) of the particle, 
d\{p) (d r (p)) is the operator of creation (destruction) of antiparticle. Sym- 
bol "f" indicates the Hermitian conjugate. Operator vf)(x) satisfies the Dirac 
equation and evolves according to the Dirac representation, satisfying the ex- 
pression 

i^P- = mx),H ]. (2) 

H in eq. (2) is a free Hamiltonian. Note that we will write 4 - vectors in 
the form A = (A,iA ). In this case the following equalities are vàlid A 2 = 
A 2 + A'l = A 2 — Al and px = p^x^ = px — Pqx ; p = Vp 2 + m 2 , m - is mass 
of the quantum of field. We use units with h, and the speed of light c taken to 
be unity (h — c— 1). Spinors u r (p), u r (—p) correspond to the partides with 
helicity r = ±1. 

We determine the operator of electromagnètic field as follows 



4 +, W=p^/^Çe; WaA ( í; )e- (i k, 

4- , w = ^/^Çe; Wa i( í; )e-- d k, 



(3) 



where k = |k|, a\{k) (a\(k)) is the operator of destruction (creation) of 7- 
quantum. Operators of the electromagnètic field A^x) evolve according to the 
Dirac representation as well. Polarization vector e^(k) is determined by the 
relations 



e (k) = (e A (k),0), e 4 (A;) = (0,i), A = 1,2,3 (polarization index). (4) 

e 1 (k), e 2 (k) are unit vectors orthogonal to each other and to the momentum 
of 7 - quantum k, e 3 (k) is unit vector directed along vector k. 
On the basis of the requirement of gauge invariance of the model the Hamil- 
tonian of interaction must be determined in the conventional form [59] (see 
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also [60]) 



#/(í) = -ie I N^(x)^(x))A^x)d^. (5) 

7 M - Hermitian 4x4 matrices (7^ + 7^ = 25^, 7+ = 7^), if> = ^74, A 7 " is 
the symbol of the normal ordering of operators, e is the charge of the electron 
so that the fine structure constant is: a = e 2 /4n ~ 1/137. 



3 Liouville formalism, "subdynamics" 



Now let me briefly examine the Liouville formalism (see for example [25], [27]). 
The time evolution of the density matrix p is determined by the Liouville-von 
Neumann equation 

Liouville operator has the form 

L = Hxl-lxH, (7) 

here symbol "x" denotes the operation (AxB)p=ApB. In accordance with 
formula (7), L can be written down in the sum of free part L that depends 
on the free Hamiltonian H and interaction part Lj that depends on Hf. L = 
Lq + Lj. Let state |a) be the eigenstate of the free Hamiltonian Ho\a)=E Q \a) 
with the energy E a . Then dyad of the states |o;)(/5| is the eigenstate of operator 
L : L \a){P\ = (E a — Ep)\a)(f3\ or L \u)) = w v \v)), where the designations 
\v)) = \a)((3\ and w v = E a — Ep were used. v is the correlation index: v — 
if a — (5 - diagonal case (vacuum of correlation) and v 7^ in the remaining 
off-diagonal case (the details of the theory of correlations can be found in 
works [2], [10], [25]). In the Liouville space for the dyadic operators we have the 
relations: inner product defined by (where symbol Tr denotes the calculation 
of the trace) 

((A\B)) = Tr(A^B), (8) 
the matrix elements are given by 

((a(3\A)) = (a\A\(3), (9) 
the biorthogonality and bicompleteness relations have the form 

({a/3'\a/3)) = (WW E l«0»<M = l - (10) 

It was shown that the description of the irreversible processes at the micro- 
scopic level is possible if eigenvalues of Liouvillian are generally complex. 
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Thus, for the Liouville operator L the eigenvalues problem is formulated for 
the right-eigenstates \^)) and for the left-eigenstates ((^J | [25] 

L|*p> = ZJ|*p>, ((fr.\L = ((^\ZJ. (11) 

Since L is Hermitian the eigenstates are outside the Hilbert space [27]. In this 
case the corresponding eigenstates have no Hilbert norm [11], [25]. For |^)) 
and (fèjl we have the following biorthogonality and bicompleteness relations 

= M«'. E = i- (12) 

Index j is a degeneracy index since one type of correlation v can correspond 
to different states. 

It is shown in work [25] that the eigenstates of L can be written in terms of 
kinetic operators C u and D v [61], [62]. Operator C u creates correlations other 
than the v correlations, D u is destruction operator. The use of the kinetic 
operators allows to write down the expressions for the eigenstates of Liouville 
operator in the following form 

= TO 1/2 *cl«i>>> = m^ D {N^I\ (13) 

where 

<^> u c = P v + C, <^> U D = P U + D v (14) 

and Nj - is a normalization constant. In the general case the operators P v 
satisfy the following condition [11] 

p u = J2 K»«^l, m<>)) = M»'- (is) 

i 

Similarly for P v and ({vj\ we have 

^ = £K»«^I, (W» = M«'- (is) 

3 

The determination of the states \uj)), {{vj\ can be found in work [25]. Substi- 
tuting (13) in (11) and multiplying P v from left on both sides, we obtain [25] 

<r c \v$) = Z»\u»)) : (17) 

where 

9 V C = P V L$ V C = L P V + P V L^ V C P V . (18) 

9q is the collision operator connected with the kinetic operator C . This is 
non- Hermitian dissipative operator, which plays a main role in nonequilibrium 
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dynamics. As was shown in ref. [11] the case v — leads (P c to the collision 
operator in the Pauli master equation for weakly coupled systems. 
Analogously it is possible to obtain equation for operator 6 V D , which is con- 
nected with the destruction kinetic operator D v . 

((v»K = ((v»\Z», (19) 

where 

6 V D = L Q P V + P V <5> V D L I P V . (20) 

Comparing eqs. (11), (17), (19) we can see that \uj)) and ((vj \ are eigenstates 
of collision operator 6qí D \ with the same eigenvalues Z v - as L. 
Thus, determination of the eigenvalues problem for the Liouville operator L 
outside the Hilbert space leads to the connection of quantum mechanics with 
kinetic dynamics. 

Operators § V D satisfy so-called "nonlinear Lippmann-Schwinger equation" . 
For the we have [25] 

= P» + Y P» — \L I & v r -& r L I ® , r \P v , (21) 

where the time ordering is introduced. For the determination of the sign of the 
infinitesimals it is necessary to determine the degree of correlation d^y 
This was defined as the minimum number of interactions Li by which a given 
state can reach the vacuum of correlation. It is assumed that the directions to 
the higher degrees of correlation are oriented in the future, and the directions 
to the lowest degrees of correlation are oriented in the past. This leads to the 
relations [25], [63] 

= +e if >d u {t> 0); = -e if d^ < d v (t < 0). (22) 

For the Q V D we have the equation 

= p» + P^ D U - ^ D Lj^ D ] £ P"— (23) 

Eqs. (21), (23) determine the kinetic operators of creation C and destruction 
D v as follows 

C " = E P * ^hï e - %L^ C ]P», 

rt» w w ^ (24) 

D v = p v w n L I - ^ d lm y P» . 



The spectral representation of the Liouville operator can be written down in 
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the form 

£ = £^K»«% (25) 
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In the Brussels-Austin group approach " subdynamics" is called the construc- 
tion of a complete set of spectral projectors IF [8], [47], [51] 

n v = EI*i»«*il- ( 26 ) 

3 

The projectors IF satisfy the following relations 

ïl u L = LU u , (commutativity); y^IF = 1, (completeness); 

(27) 

ITir = IT<U,, (orthogonality); IT = (IF)*, (star-Hermiticity), 

where the action "*" corresponds to the "star" conjugation, which is Hermi- 
tian conjugation plus the change — > e Ufl [10] [11]. Oprerator IF can be 
represented in the following form [25] 

IF = {P v + C V )A V {P V + D V ), (28) 

where A" is the star-Hermitian operator 

A v = P\P V + D V C V Y 1 P V . (29) 

Taking (27) it is possible to write down the density matrix p as follows 

p = £lFp = ]>>, (30) 

V V 

where p u = IF p. Projectors W can be associated with the introduction of 
the concept of "subdynamics" because the components p v satisfy separate 
equations. In the framework of the "subdynamics" approach the time evolution 
of the density matrix has the form (see, for example, work [11]) 

fT(t) = n>(í) = exp(-iLt)Wp(0) = 
= (P" + C v )e- ig c*A v (P v + ir)p(0). [ } 



4 Time evolution of the density matrix 



Multiplying IF and P v from the left on both sides of eq. (6) we obtain 

^7 = P v Lp". (32) 
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Let operators Q v determine subspace ortogonal P v . The states belonging to 
subspace Q v have a degree of correlation differing from those which have the 
states belonging to subspace P v 

P»Q" = Q"P" = 0, (33) 



P" + Q" = 1. (34) 

Using oprerators P v and Q v we can rewrite eqs. (24): C v = Q V C V P V ', D" = 
P V D V Q V . It is easy to see that operators C v describe transitions from P v 
correlation subspace to the Q u correlation subspace and operators D v describe 
transitions from the correlation subspaces other than v to the v subspace. The 
expressions (28), (32) and (34) result into [25] 

i^py = w . (35) 

The case u=0 leads (35) to the kinetic Pauli master equation for P v p v - com- 
ponent. Eq. (35) describes the time irreversible evolution of the unstable state. 
Our great interest is to investigate eq. (35) for the system of interacting rela- 
tivistic Fermi - Dirac particle and electromagnètic field. For this purpose we 
will obtain the Liouville - von Neumann equation for P v p v - component in the 
Dirac representation. In the Dirac representation for eq. (35) we have 

i^P u p u (t) = P»L I {t)C v P v [P» P v {t)\. (36) 



Determining operator 



for eq. (36) we get 



r (t) = P v L I {t)C v P v (37) 



^P7(í) = r(í)P7(í), (38) 

where the previous designations of operators are preserved. 

The general solution of eq. (38) can be found after examining the equivalent 

integral equation. The solution of eq. (38) we will search for the component 

py(t) 

P u p u (t) = P u p u (t ) + H) / dt 1 ·& v {t 1 )P v ( f{t 1 ), (39) 

to 

where P u p u (to) is the component into initial t — to moment of the time. 
Substituting in the right side of the expression (39) instead of P v p v (ti) the 
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sum 

íi 



py(ío) + H) / dt 2 {r{t 2 )p v P v {t 2 ) (40) 



to 

and consecutively continuing this procedure we find 

í t íi 

p v p v (t)= (i + (-i) | díií?"(íi) + H) 2 1 1 dhdhritjrih) + 

ío ío ío 

í íl tn-1 (4]^ 

... + (-i) n f I ... I dt 1 dt 2 ...dt n r(t 1 )r(t 2 ) ...r(t n ) + ...)x 

ío to to 

X PV(ío). 

The obtained expression can be written down in the form 

P"p"(í) = ÍT(t, ío)PV(ío), (42) 

where 

íí"(í, to) = EH)"/ / - / «2 ... dt n r(h)r(t 2 ) ... r(t n ). (43) 

n=0 to t ío 

Operator Q^í, t ) i n expression (42) plays the role of the evolution operator. 
The non-Hermitian operator Q u (t, t ) determines the time irreversible evolu- 
tion of the density matrix - the time irreversible evolution of the relativistic 
unstable state. 



5 Complex spectral representation of Hamiltonian 



Let me examine the eigenvalues problem for the Hamiltonian H = H + Hj. As 
for the Liouville operator the problem will be formulated outside the Hilbert 
space. In this case as earlier we must distinguish equations for the right- 
eigenstates |</? 7 ) and for the left-eigenstates (</5 7 | of Hamiltonian, where 7 is 
the index of the state 

(if + #j)K) = £> 7 ), (^\(H + Hj) = (^\Z^ (44) 

eigenvalue Z 7 is the complex number. Since H is Hermitian the corresponding 
eigenstates |</? 7 ), (</? 7 | have no Hilbert norm. 

< </? 7 | </? 7 >=< </? 7 | </3 7 >= 0. (45) 
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The "usual" norms of the states |</? 7 ), (</? T | disappear as required to preserve 
the Hermiticity of H [10] (the details of the complex eigenvalues problem can 
be found in works [2], [26]). Let me write down the Hamiltonian of interaction 
in the form Hj = gV, determining explicitly coupling constant g. The value of 
coupling constant depends on the model of interaction and will be determined 
later. Solutions of the eqs. (44) can be found after presenting vàlues |</? 7 ), (y5 7 |, 
in the perturbation series 

oo oo oo 

M = £<? n K n) >, (*M = £í? n <^ n) l, z, = 5>»z<">, (46) 

n=0 n=0 n=0 

where 

K 0) > = |7>, (^ 0) l = (7|, Z^ = E T (47) 
As it was shown in ref. [10] from relations (46), (47) we can obtain 

n— 1 

= (l\V\^- 1] ) - - £ ^(tI^), (48) 
i=i 



w\- ((/^l^) - (49) 
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where in accordance with Brussels - Austin group approach the time ordering 
is introduced. The sign of infinitesimal depends on the direction of the 
processes: the transition 7 — > (3 we will associate with = e > 0. 
Define I7) as a bare state, which corresponds to relativistic Fermi-Dirac parti- 
cle and \f3) as a state consisting of the bare Fermi-Dirac particle and photon: 
I7) = |p, r), = |p, r';k, A), where |p, r) refers to a one - particle state, 
|p',r';k, A) is a two - partides state, p (p'), r (r') - momentum and helicity 
of the particle and k, À - momentum and polarization index of photon. In the 
model, states |p, r) (|p', r')), |k, À) are eigenstates of the free Hamiltonian Hq\ 
Ho\p,r) = E p \p,r), íío|k, A) = c<Jk|k, A) with E p = TpTn? and uj^ = |k| 
(m - Fermi-Dirac particle's mass). Hamiltonian Hj is determined by the ex- 
pression (5), g = e is the charge of the electron. Substituting the expressions 
for ip(x) (1), A^x) (3) in (48), multiplying by e n and summing with respect 
to n, we obtain 

^ m í v-^ rfp'rfk' , , .. 

^ = ^ - « ^5372 y £ (w^yv5 *<p - p + " > - 
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The expression for (p', r'\ k', X'\(p p , r ) can be obtained from relation (49). In 
our case we get 

(p , r ; k , a w v , r ) = «— j E ( ^ p ,, 2av)1/2 x 
x ^VW) 

where -Ep^k' = E p > + uv. Substituting the expression (51) into (50) for Z p , r 
we obtain 

^{k') 7 ^{p)lixU r '{p')(íy '(k')u r '(p')'j u u r (p) 
E p ' t k' — Z pr — ie 

Using the formal expression — > p-^ =|= m8{w) rewrite (52) in the form 

Z Pt r — E pr — ï7 Pi r, (53) 

where 



(54) 



E p ' : k> — Z pr 

is the renormalized energy (p stands for the principal part) and 

m 2 r „ dp'dk' 
Tp, = -^i5,W^ í(p - p + k) X (55) 

X e^V(pbXV)eí'(*VVhX(ï>) 1 5(£p.,k· - Z p ,r). 
The expression (52) leads to the relation 



(56) 



|<(fr)^(pbxV)l 2 _ n 

X ~~ ~~ ; — U. 

From eq. (56) we can obtain the connection 

E n l b-/ — Z n r — %S E-r.1 u/ — E-r, — %E ' ^ ^ 
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where 0(e) determines the terms of higher orders on e. Limiting the expres- 
sion (55) by order e 2 we have 



m 2 r ^ dp'dk' . , u , n ,_, . 

T~ - ^ / E «(P " P + « )í(^ - B P ) x (58> 

where for the photon À' = 1, 2. Since the expressions (56) (57) contain two 
complex terms Í7 Pi r and ie, which determine a pole in the lower half plane and 
a pole in the upper half plane, the operation of integration of the expressions, 
which contain the vàlues of the form ■= h= is accepted to determine as 

follows [10], [11], [48]: we first have to evaluate the integration on the upper 
half-plane (C + ) and then the limit of z — > — ï7p, r must be taken. For example 
the integration over E p/ with a test function g(E p >) can be presented as follows 



z 



hm (ídE p 9 -^l 

^-iip,r\J E v ' + U]ç> — E v — zJzec+ 

oc (59) 

, í iEp m . . 

{ P (£V + "V -E v - z)t ilt r 

This special feature will be used below for the determination of the expression 
for the density matrix. 



6 Expression for the density matrix 



Determination of the expression (42) we will carry out for the diagonal (^=0) 
matrix element of the form ((77|p°(í))) = P° 7 (0- Let initial moment of time 
be zero. In this case the expression (42) can be represented as follows 

t 

p 77 (í)=p 77 (o) + (-<)/díi«Ty|^(*i)l««»pL(o) + 



t tl 

+ H) 2 / /rfí 1 rfí 2 ((77l^ (íi)l7V))((7 , 7V°(Í2)|«a))pL(0) + 



o o 

í tl Í2 



+ H) 3 / / /rfí 1 «3((77l^ (íi)l7V))((7'7V (í 2 )l7V)) x 

x ((7 / Yl·°(í3°)l««))pL(o) + 



(60) 



where summation (integration) over all infernal indices 7', 7",.., a is implied (It 
is necessary to note that for the simplification of the expressions the normal- 
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izing volume is implied, but it is not written.). Eq. (21) leads to the following 
approximation for the collision operator Q v c 



9 U C ^L P U + P V L 1 C V P\ 



where 



-1 



u^u wn-W- xe^ 



L T P V . 



(61) 



(62) 



From the expressions (37), (62) it follows the expression for the operator i?"(í) 



r (í) = J2 P u Li(t)p» 



-i 



— w u — ie, 



(63) 



Then for the i?°(í) taking into account the condition > d we obtain 

-1 



0°(t) = J2 P%(í)P A 



— ie 



L!{t)P . 



Using result (64) we examine the second term of expression (60) 



(64) 



-i) / ^i((77l^°(íi)l«a))pL(0) 



{-i) / dh 



a 
t 



(7|^ J |a>(a|^ / | 7 ) H#/|7>(7l#/|a> 



w a ~i - ie 



+ 



w^ a - ie 



PL(0) - (65) 



- (-o / dt x J 



< 7 |tf,|p><p|#,|7> (p|///|7)(7l^/|p> 



+ 



/4(o). 



o P 

Symbol / in (65) indicates summation over discrete and integration over 



continuous variables. Selecting the state (7) in the form |p, r) and taking into 
account the expressions (5), (57), (59) we obtain 



(7|^/|a)(a|^|7) {^H^^H^a) 



w a ^ - ie 



+ 



w^ a - ie 



PL(0) = 2 7p , r e 2 x 



x / v dp'dk' < p ' r l^ p/ ' r/; X'WS-' 



(66) 



(3 



P ',k' 



E r 



E r 



X Pp',r';k',A' p',r';k',A' (0) ) 



where U = J N(ifj(x) , y l·l ifj(x))A l·l (x)d^. and 7 Pjr is determined (55). The desig- 



nation (-Ep',k' — E p 



-*7p, 



corresponds to the integration, which first of all 



is carried out in the lower half complex plane (C ) and, after that, the limit 
of z — > +ï7p, r is taken. 
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For convenience of the further consideration let me introduce the new desig- 
nations. We define the function 



F TÍ = 2 7r e 2 



(r\U\i)(i\U\r) 



( Ei - E T - z)ti^(Ei - E T - z) 



(67) 



>+Í7t 



where Greek and Roman indices \r), \i) correspond to the one-particle and two- 
particles states, respectively. For the function T ri it is possible to determine 
the rules 



r„p°(0)=2 7r e 2 j 



(r\U\i)(i\U\r) 



r„r /3í p°(0)=2 7r e 2 J 



(Eí-Et-z^ÍEí-Et-z)^ 
(r\U\t)(t\U\r} 



P°(0), 



x 2 7/3 e 2 



(Ei -E T - z)t ilT (Ei E T — z) +ilT 
((3\U\i)(W\P) 



x 



(Ei-Ep-z^M-Ep-z) 



-p°(o). 



(68) 



The rules (68) lead to the following expression 



H) 



+ 



w a "< - ie w^ a - ie 

a 

r p , r p',r';k',A'Pp' jr ' ; k',A' p',r';k',A' (0) ■ 



pL(o) = 



(69) 



Substituting (1), (3), (5), (57), (59) in the second term of expression (65) we 
obtain 



H) 



- ie 
p 

27 P ,rP p , r p , r (0). 



+ 



w^p - ie 



P? 7 (0) = 



(70) 



In expression (70) integration (summation) over the state |p, r) is not carried 
out. Finally we obtain the result 



H) / rfíi(( 77 |^°(íi)|aa))pL(0) = (r p , r p ,,r> ; v,x x ^ 

X Pp',r';k',A' p',r';k',A'(0) ~~ ^1p,rPp,r p,r(0)K 
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Analogously for the the third and fourth contributions to the expression (60) 
we obtain 

t íi 

H) 2 / /rfí 1 rfí 2 ((77l^ (íi)l7V))((7 , 7V°(t 2 )|aa))pL(0) = 
o o 

= ((27p,r) 2 Pp,r p,r(°) ~ 27p,rF p , r P ',r';k',A'P p /, r / ;k ',A' p',r';k',A' (°) + (72) 
+ Tp jr p ',r';k',A'r p ", r " p',r';k',A'Pp", r "; p",r"(0) _ 

/t 2 
J2dp"T p „y, p',r';k',A'Pp', r ';k',A' p',r';k',A'(°)) > 



t íi í 2 

H) 3 / / /rft 1 « 3 ((77|^ (íi)|7V))((7 / 7V (Í2)|7 ,, 7")) x 
x ({l"l\l\h)\aa))plM = (-(27p, r ) 3 p P ,, p>r (0) + 

+ (27p,r) 2 r p , r p',r';k',A'Pp', r ';k',A' p',r';k',A'(°) + 27p,rF p , r p',r';k',A' x 
X / rf P"r p '',r" p',r';k',A'Pp', r ' ; k',A' p',r';k',A'(°) + ^p,r p',r';k',A' X 

X Tpiiyi p\ r ';k',X'^p",r" p"',r"';k"',A"'Pp"' ) T·"';k"',A"' p"',r"';k"',A"' ~~ 

— r p , r p',r';k',A'r p " jr " p',r';k',A'27p",r" x (73) 

X P p »,r" P ",r"(0) _ T Pjr p ',r';k',A' / X ^p'TpV" p',r';k',A' x 
x r p /// ;7 ./// p',r';k',A'^p'",r'" P "',r"'(0) + ^P,r P ',r';k',A' / X '"fc" X 

/ XI ^P W r p ''',r'" p',r';k',A'P p ' ir ' ;k ',A' p',r';k',A' (°) ~ 



X r p // pfy.^^x' 

o f3 

- 27 P)r r Pjr p ', r ';k',A'r p ", r » p',r';k',A'P p ", r » p»,r"( )) gjJ 

where the procedures of integration and summing are achieved on all contin- 
uous and discrete repeating indices besides the indices p, r which correspond 
to the state I7). 

Studies of the expression (60) lead to the genealogical connections, where 
each of the foregoing contribution gives birth to the following contribution 
which determines the sequential term of the sum. For example, contribu- 
tion — 27 7 p° 7 (0), which determines second term in the expression (60), is 
the ancestor of the contributions (27 7 ) 2 p 77 (0), — 27 7 r 7i p° i (0), -(27 7 ) 3 p 77 (0), 
(27 7 ) 2 r>°(0), ...,(27 7 )> 77 (0), -(2 77 )™- 1 r 7 ^(0), ... . Analogously, contri- 
bution r 7 ip°j (0), determining second term generates contributions: r^r^pJg^O), 

— r 7 j / T pi p% (0) and so on. For example, for the fifth order we have the con- 
p 

nections 
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(27 7 )V° 7 



-(2 77 ) 5 r 7íP o i -(2 77 ) 5 r 7í r /3í p° /3 (2 77 ) 5 r 7í /r^ 



•(2 77 ) 5 r 7l r /3íP o /3 



(2 77 ) 5 r 7l 2 7/3 r /3íP o /3 _(2 77 )5r 7i r /3i r^ P 2fc 



where integration and summing are achieved on all repeating indices besides 
7. Further analysis of the terms of the expression (60), leads to the extremely 
great variety of contributions of higher order on e. I limit my analysis by the 
contributions determining the structure of the density matrix in the approxi- 
mate form 



5° 
P,r P,r 



(*) 



-2 7p , r í 



P° P , r p, r (0) + (1 - e- 2 ^)T p , r pV , ;k , y x 



X Pp',r';L·',X' p',r';k',A'(0)- 



(74) 



The expression (74) is determined so that the function T p r p>y,k',\' does not 
contain the value 27 Pjr . Expression (74) follows from equation (35) and cor- 
responds to the kinetic, irreversible evolution of the unstable electromagnètic 
system in the time to the equilibrium state. Thus, the ordering in the time 
leads to the complex eigenvalues. Such complex eigenvalues make it possible to 
describe the relaxation process in other words the irreversible process without 
appearance of the other spontaneous, unstable states. 



7 Numerical calculation 



I examine the first term of the expression (74) which is the probability of 
finding Fermi - Dirac particle with momentum p and helicity r depending on 
the time. We will assume that the particle is not polarized. The averaging over 
r leads to the following result 

Ve- 2 ^Vp,, p>r (0) = e-^^Wplr P ,M = Pl M (75) 
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where 

m 2 dp'dk' , , 

7 — ^ - e V/Çw^ í(p - p+k)x 

x 5{E p ,^, - E p )^Tr(e x '(k') ■ yA(p')e x ' (k') • 7 A(p)) ( 76 ) 
p + im 

and A(p) = ~2Ï^T" 

Using the expression (28) we represent the density matrix p pr pr (0) in the 
form 

P° P , r p, r (0) = P P ,r p, r (0) - (D°C°p(0)) p , r p>r . (77) 
From the expression (62) and relation = (C^)* [25] we can find 

Plr P>r (0) - Pp, P ,(0) - e 2 7^ / E Í(P' " P + k') x 

where c.c. means the complex conjugate. Density matrix p PjT p , r (0) has the 
form [59] 

Pp, P ,r(0)=« r (pK(p). (79) 

We determine the expression for p p (t) being limited to term of lowerst order 
on e. In this case we have 

pl p (t) » e-^-^Vp P (0) (80) 

with 

ft.p(o) = ^E« r (pK(p) = ^A(p)- (si) 

relativistic density matrix of the Fermi-Dirac particle. Let estimate the den- 
sity matrix summing up the diagonal elements of the expression (80). This 
procedure results into 

(P P P (t))dia 9 ~ e- 2 ^^\p p p (0)) diag , where (p p p (0)) diag = 1. (82) 

Since 7p, T ·=±i depends on the momentum p we examine the special case, when 
the angle d p of vector p (in spherical coordinates) is zero. Summation over A' 
and integration over 5 - functions give the result 

Oí 

7| P |,r=±i = — (/1 + /2 + /3 + /4), (83) 
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where 

= _^PP f cos 2 (fl k /)tft2 k , (M) 

2E P J (\p\ 2 + ou 2 + m 2 - 2\p\üü cos(iï kf )y/ 2 + ou - |p|cos(tf k ')' 1 ' 



cj 2 |p| r cos(^k')^k' 



2E P „ 


1 ( p 2 + (jj 2 + m 2 ■ 


- 2 p cucos^k')) 1 / 2 + w — P cos(i?k') ' 




(| P | 2 + u 2 ^ 


- m 2 — 2 p cjcos(^ k ')) 1 ^ 2 ^k' 


-lí 


( p 2 + u; 2 + m 2 — 


2 p \lu cos(^k')) 1 ^ 2 + 00 ~~ P cos(í9k') ' 


um 2 


/ 




2E P 


i ( p 2 + uj 2 + m 2 


— 2 p w cos(-^k')) 1 ^ 2 + u ~ P cos(i?k') 



35) 



u/ ; ^|p| rw -riu — z|p|í^uus^í/k';; ' u,"k' / c «x 
y 3 = 77 / , .9 , ...o ni_i /.o nm/o , . .i _/.o n > W 



(87) 

and is the energy of photon (without the Doppler effect. In our case the 
Doppler effect is not significant). The calculation of the expression (82) was 
accomplished numerically with the use of program Mathematica and the fol- 
lowing approximation 

n 12 2 2 «i i /n \\l/2 |p| 2 + Cc> 2 + m 2 - Iplcj COs(-# k ') / ooN 

p 2 + ^ 2 + m 2 -2 p ^cos tf k ' 1/2 ~ 7i i^— 2 ' 2W2 88 

(|p| 2 + w 2 + m 2 ) 1 /^ 

The calculated results are represented in figs.1-3, where the process like e — > 
e + 7 (bremsstrahlung of electron) is examined. The calculations are executed 
for the different vàlues of momentum |p| of electron. It is seen, the time evo- 
lution of the density matrix depends on the value of momentum of particle: 
the density matrix is decreased with an increasing of momentum |p|. 
It is necessary to note that these calculations do not determine the physical 
bremsstrahlung. First of all because the used model of interaction is deter- 
mined for the bare not dressed partides. Furthermore it is known that electron 
in the bare, free state cannot either absorb or radiate the photons. Strictly 
speaking, the energy uj of the radiated photon in the expressions (84) - (87) 
must be zero. Nevertheless, in the model it is assumed that electron interacts 
with the external electromagnètic field. Therefore, the energy uj of the radiated 
photon was forced different from zero. The model makes it possible to develop 
the procedure for the description of the realizable irreversible processes. In the 
work [64] in the framework of Prigogine's principies the weak interaction like 
n ± - meson decay is investigated. 



8 Concluding remarks 



Let me briefly summarize the results. The time irreversible evolution of the 
relativistic, unstable electromagnètic system is investigated in the framework 
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of Prigogine's principies of description of nonequilibrium states on the ba- 
sis of unified formulation of quantum and kinetic dynamics. As a result the 
expression for the density matrix determining irreversible evolution of the rel- 
ativistic unstable system in the time was obtained. Although I do not examine 
the question of the determination of observed physical process, the approach 
makes it possible to define the expression, which can be initial for the further 
construction of the irreversible relativistic model of time evolution of the phys- 
ical relativistic unstable systems. It is interesting to investigate the possibility 
of applying the developed procedure for the time irreversible description of 
the observed physical processes such as relaxation of the unstable states of 
àtoms and atòmic nuclei, bremsstrahlung, particle decay. All these problems 
are very debatable and require further consideration. 
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L·t 



2 -10 5 4-10 5 6 -10 5 8-10 5 

Fig. 1. Density matrix (p p p {t)) diag : m = 0.51 MeV, |p| = MeV, lü = 12.8 eV, 
t(sec.) 




Fig. 2. Density matrix (p° p (t)) d iag ■ m = 0.51 MeV, |p| = 0.001 MeV, uj = 12.8 eV, 
t(sec.) 




t 



2-10 4-10 6-10 8-10 



Fig. 3. Density matrix (p° p (í))di og : m = 0.51 MeV, |p| = 0.01 MeV, w = 12.8 eV, 
t(sec.) 
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